Thermal transport of molecular junctions in the pair tunneling regime 
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Charge and heat transport through a single molecule tunnel-coupled to external normal elec- 
trodes have been studied. The molecule with sufficiently strong interaction between electrons and 
vibrational internal degrees of freedom can be characterized by the negative effective charging en- 
ergy U< 0. Such a molecule has been considered and modeled by the Anderson Hamiltonian. The 
electrical conductance, thermopower and thermal conductance of the system have been calculated 
as a function of gate voltage in the weak coupling limit within the rate equation approach. In the 
linear regime the analytic formulae for the transport coefficients in the pair dominated tunneling are 
presented. The effects found in the nonlinear transport include inter alia the rectification of the heat 
current. The sense of forward (reverse) direction, however, depends on the tuning parameter and 
can be controlled by the gate voltage. We also discuss the quantitation of the thermal conductance 
and the departures from the Wiedemann-Franz law. 

PACS numbers: PACS numbers: 73.23.-b; 73.63.Kv: 73.23.Hk 



o 



> 

ON 

in 
o 
m 

o 



X 



I. INTRODUCTION 

The appearance of the effective attractive interactions 
between electrons in metals^ is responsible for the insta- 
bility of the Fermi surface and the resulting phenomenon 
of superconductivity; one of the most spectacular quan- 
tum effects on a macroscopic scale. The strong local elec- 
tron - phonon interaction may induce formation of po- 
laron or bipolaron quasiparticles which are at the heart 
of bipolaronic theory of superconductivity^i^. The so- 
called negative U centers^ are not only a source of super- 
conducting instability but also play an important role in 
physics and chemistry of materials^ (see, however—). 

The tendency towards bipolaron formation is 
enhanced^ in the confined structures like quantum 
dots^. This assertion has been corroborated by means of 
quantum Monte Carlo calculations in the strong coupling 
regime for various dimensionalities of nano-structuresS. 
Negative charging energy of the small systems like 
quantum dots is a relatively new concept*i2rJ^. Very 
small system is characterized by the small capacitance 
C which makes the charging energy e^/2C large. To 
overcome it sufficiently strong polaronic shifts resulting 
from strong coupling of vibrational and charge degrees 
of freedom in a molecule are needed. The charge fluc- 
tuations in a device with metallic and superconducting 
grains may also over-screen the Kondo repulsiorti^ in the 
normal grain. 

The purpose of this paper is to extend the recent 
worki^ in which charge transport through the molecule^ 
characterized by the negative charging energy has been 
studied. Modeling the system using Anderson like Hamil- 
tonian with a negative value of the charging energy U and 
considering weak coupling between the molecule and the 
electrodes the authors have systematically accounted for 
the processes contributing to charge transport and calcu- 
lated its differential conductance. The study of the neg- 
ative U molecule has later been extended to the stronger 
coupling and the regime of the charge Kondo effeclii^. 



Here we consider the same model and study the charge 
and heat transport via negative U molecule in the weak 
coupling regime when the pair tunneling processes domi- 
nate the transport. The analytical formulas for the ther- 
mopower S and the thermal conductivity k in the linear 
regime of the low bias voltage V = {^l ~ Ii-r)/^ and the 
temperature difference ST = Tl — Tr have been found. 
t^L{R) denotes the chemical potential and Ti^[r) temper- 
ature of the left (right) external metallic bulk electrode, 
e denotes the positive electron charge. As in the previ- 
ous work we assume that the molecule is characterized 
by single, doubly degenerated electronic level. In the 
nonlinear regime with large temperature differences 6T 
and asymmetric coupling to the external leads the system 
shows heat rectification property. The thermal current 
and conductance in the forward direction (say, ST > 0) 
differ from those in the reverse direction (ST < 0). 



The low temperature thermal conductance k in the 
linear regime is proportional to T, but the heat quan- 
tum, defined as k/T takes on a non-universal value. 
However, the Wiedemann-Franz law is obeyed in the 
limit r — > 0. At elevated temperatures we observe 
marked departures from the Fermi liquid behavior with 
the ratio L = k/{GT) exceeding the Lorenz number 



in 



The organization of the rest of the paper is as fol- 
lows. In the next section we briefly recall the model 
and approach. In section HI we present the transport 
coefficients calculated in the linear regime. In section 
IV the nonlinear transport coefficients are defined and 
calculated. In section V we discuss the validity of the 
Wiedemann-Franz law and quantization of the thermal 
conductance in the context of pair tunneling. We end up 
with summary and conclusions. 
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II. THE MODEL AND APPROACH 

The aim of this section is to recall the main results ob- 
tained by Koch et a/.^S and to establish the notation. 
One starts with a general system consisting of exter- 
nal leads and a central molecule coupled by tunneling 
processes. The electron subsystem on a molecule inter- 
acts with bosonic degrees of freedom (phonons), which 
are eliminated by performing the Firsov-Lang canonical 
transformation. We assume polaronic shifts large enough 
to make the effective charging energy Ue/f = U — 2X^haj 
negative. We denote it simply U. The motional nar- 
rowing of hopping parameters as well as shifts of the dot 
energy level are absorbed in the definition of the corre- 
sponding parameters in the Hamiltonian 



H — Hmol + Hi 



leads 



(1) 



where Hmoi = ^dnd + Und-\ndi describes the interacting 
electrons on the dot. Here nd — ^da — J^a^a'^'^ 
is the number operator, and d\{d„) denotes the cre- 
ation (annihilation) operator of a spin a electron on 
the molecule. The external leads (left-L and right-R) 
are characterized by the electron energy spectrum et. 
The corresponding term in the Hamiltonian is given 

by Hieads = EA=L,_R,Ek,<TCA,k4kaCAk<T, with k = 

£a k — Ma and jix = ji — eV\ denoting the chemical po- 
tential of the A electrode subject to the bias voltage V\. 
In the following we assume the equilibrium value of the 
chemical potential n = and e^k = Cflk — Ck. 

The coupling between the molecule and leads is gov- 
erned by the Hamiltonian 



H, 



\=L,B. k.cr 



Aclkc 



h.c). 



(2) 



To correctly account for all low energy processes which 
contribute to the transport in the limit of negative U, it 
is advisable to eliminate Hi by means of the SchriefFer- 
Wolff^^ transformation. This is valid in the limits I^aI ^ 
\U + ed\ and \t\\ ledT''. One gets^ the effective low 
energy Hamiltonian 



H — Hmol + Hieads 



Hn, 



(3) 



in which the direct and exchange Hdir,ex interactions be- 
tween the dot and the leads read 



1 V- 



1 



^Ak — £d 



'^Akcr'^A'k'cr 



(4) 



AA'kk'cr 

M{exk){di^dacli^„Cyu' -a - 4ko-CA'k'o-n-d5-) + H.C. 



where M(e) = [e — £d\~^ ~ [s — £d — U]~^ . These terms 
are most important in the study of transport through the 
quantum dots with positive charging energy, but they 
do also play a role in the present case of 'negative U 



quantum dot'^^ and describe single electron co-tunneling 
processes. 

The pair terms readiS 

Hpair^ ^ txtl,M{ex\^)d^d^c\,^^,^c\^^ + h.c.. (5) 

AA'kk' 



For negative values of U , the pair tunneling terms play 
a main role because the state with two electrons on the 
dot is its lowest energy state and the double occupancy 
of the dot is favorable. The two electrons may tunnel 
onto the dot from a single lead or from two different 
leads. The single occupation of the dot is not favored 
for the low voltage bias eV <^ \U\ and the sequential 
single particle processes are exponentially suppressed. As 
a result the single particle events can contribute to the 
transport through negative U molecule via higher order 
processes (co-tunneling) . 

The approach is based on the rate equations for the 
occupation probability P{n), where n denotes a state of 
the molecule. Since for negative U the single occupation 
of the molecule is never favorable, one finds either one or 
two electrons on a molecule. The energies of these two 
states are equal if 2ed + U ~ 0. Due to the normaliza- 
tion condition P{0) + P{2) = 1 the rate equations reduce 
to dP{2)/dt = P(0)VKo_^2 - F(2)PF2^o- In a stationary 
state this gives P(0) — ^^2^0/ (W'o-»2 + W2^q), where, 
in the notation of the paper [Ifll , W^i-s-/ is the total tran- 
sition rate from the initial state i to the final state /. 

To calculate transition rates we use the Fermi's golden 
rulei^ and consider all processes in which electrons or 
electron pairs tunnel from the electrodes A and A' to the 
dot or vice versa. With a constant temperature across 
the system they have been obtained by Koch et al. and 
for the process —J' 2 readi^ 



W, 



AA' 



TaFa' 



0^2 - —r- I deM\e)f{e~eVx)f{2ed + U-e-eVy) 

(6) 

Equation ^ describes the transfer rate for tunneling of 
2 electrons onto the dot (spin up electron hops on the 
dot from the electrode A and spin down one from the 
electrode A'. In the formula © = 2TTpx\t\\'^ denotes 
the effective coupling between the dot and the lead A, 
and Pa stands for the density of the states at the Fermi 
level in the electrode A. fx{x) is the corresponding Fermi 
function. Similar expressions can be derived for all other 
rates. The rates of single electron co-tunneling from the 
left to the right electrode, without changing the occu- 
pancy of the dot i.e for the processes 0^-0 and 2-^2, 
are obtained^° from Hdir,ex- Obviously Eq. ^ is also 
valid for system in which the left lead is characterized 
by temperature Tl different from that of the right lead 
Tr. In this case, however, no general analytical solution 
is possible, except in the linear regime. It is presented in 
the next section. 
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III. THE LINEAR REGIME 

First let us consider the transition rate of an electron 
from a given state in the electrode A into the quantum dot 
as induced by the co-tunneling term in the Hamiltonian. 
Let the initial state \0)d of the dot be the vacuum |0) (no 
electrons on the dot). The initial state of electrodes rep- 
resents two Fermi surfaces with the electrons occupying 
the states up to ^LifJ'R)- We denote the initial state of 
the whole system by \i) and its energy by Eq. The pro- 
cess in which an electron with the spin a from the state 
k in the electrode A is transferred to the state k' in the 
electrode A' results in the final state |/) = c\k, 



The energy of the final state is Ej = Eo+ G'fc' ~ £.xk- In 
the single process the energy transported per unit time 
from the electrode A is ^Ak- Its contribution to the total 
energy flux through the left junction equals 



Wr 



;A,A' 



0^0 



(l-ZA'CevkOMlCAk-evk') 



/A(^Ak) (7) 



/A(£Ak) = l/(e('^^k-eV^)/fcBT^ ^ 2) is the Fermi distribu- 
tion function. 

In fact the same factor contributes to the heat flux 
in an elementary process in which e.g. a pair is hopping 
from the states Ak, A'k' onto the dot, ii X — L and X' — R 
. In this process only the member of the pair with the 
quantum numbers Lk transports the energy flux through 
the left junction. If both electrons stem from the left 
electrode then their contribution to the total flux reads 



.XX' 



27r v-^ , 



ntA'nM(e,k)p(eAk 



6'k') (8) 

/A(eAk)/A'(£A'k') S{2ed + C/ - CAk - 6'k') ■ 



kk' 



The conservation of energy expresses the fact that if ini- 
tially the energy of the system is Ei — Eq then final en- 
ergy is Ef = Eq + 2ed + U — £_xk — £,x'k' ■ In calculations of 
the current flux the extra factor of 2 appears, account- 
ing for a double charge carried in the above processiS. 

The net charge and heat currents (each of them being 
a sum of pair and co-tunneling contributions) in the left 
junction are given by 



rtot.L 
>2 



,tot,L 
O-s-2 



(9) 



W 



while M^*' 



tot 



with W^°l''^ = 2W^-^ +W'^-S: - 
rLR _ W,{^„ and completely analogous expressions for 



n—yn 



Ti/£,tot,L J we, tot 
''n— >m ^^^^"-^ yy n^rf 

In the linear regime V\ 0, ST\ — > and setting 
Ma = M + gVa we express the transition rates up to the 
linear order in V and ST, using 

h{^)~ f{^) + -Q-^eVx^—e—. (10) 



For the asymmetric couplings F/, ^ F/f the resulting for- 
mulae are long and will not be reproduced here. In the 
special case F^ = T ^ = F and for the symmetric dis- 
tribution of voltages Vl/^ — ±.V/2 and the temperature 
difference Tl/r. = T ± 5T/2 we obtain 

/e = G^V + l'^dT 
Iq = I^V + KST. (11) 
The contributions to the charge and heat flux read 



Go — 2FLFi^ — 
h 



1/2(0) 



2 sinh I3x 



K = 



(sd + uy 

2TlTr 



{Ed + uy 



+ Af2(0)- 



1 
T 



(12) 
(13) 
(14) 

(15) 



Iq is explicitly 



i sinh f3x 

In the above formulae we introduced the notation a; — 
2ed + U . X measures the distance from the degeneracy 
point. It can be tuned by the gate voltage. Note, that the 
Onsager reciprocity relation^'^ 
fulfiUed. 

It is easy to obtain the phenomenological transport 
coefficients from those currents. The linear conductance. 
Go, is defined for 5T = Q as le — GqV and is given by 
equation . It has been discussed earlier— and we shall 
not discuss it here. The thermopower S is defined as the 
voltage, which appears across the system subject to the 
temperature gradient in the absence of current fiow 



(16) 



S has been studied earlier^"^ and the analytic formula for 
it reads2i 

kn M2(0)/32a;V4sinh/3a; 



5* = -- 



Go 



(17) 



Thermal conductance k is defined as a coefficients be- 
tween the heat current Iq and the temperature gradient, 
under the condition of no charge current 1^ = 0. With 
this definition, we find 



+ 



2F,F^^ 
h 

M2(0)/32r 



fix) 



e 



8 sinh f3x 
In the above formula I2, 



S 



d i^d 



/3/2 



4 sinh Px _ 
denotes the integral 



(18) 



h = 



df 



4 
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FIG. 1: (Color online) Dependence of the heat current Iq on 
X = {2ed + U)/\U\ for (5T = ±0.003|I7| and Fl = 2,rL-rR = 1 
in units of \U\. The temperature T — 0.005|(/| and Jo — 
e\U\/h. 



IV. THERMAL RECTIFICATION IN A 
NONLINEAR TRANSPORT 

Outside the linear regime, i. e. for finite voltage bias V 
and temperature difference between the electrodes, the 



transfer rates W^^^ and W^^^ have to be calculated 
numerically. In the studies of charge conductance^ in- 
teresting rectification properties of the device have been 
predicted. For strongly asymmetric coupling of the dot 
to the electrodes e.g. Tl » the current /g is 
the asymmetric function of x in the limit of \eV\ » 
ksT . The intriguing question arises if similar rectifica- 
tion properties'^ could be achieved in the heat flow. To 
answer this we have calculated the currents /g and Iq for 
general values of V and ST. Formally both charge and 
heat currents depend on a voltage and a temperature 
difference 



h=Ie{V,5T) 

Iq = Iq{V,5T). 



(20) 
(21) 



As stated earlier the thermal conductance k is gener- 
ally defined by Iq = —k8T under the condition 7e = 0. 
The condition of no current flow defines the thermoelec- 
tric power 5" = — where V is the voltage gener- 
ated in the system by the temperature gradient. It pro- 
vides the 'selconsistent value' of the voltage to be used 
in equation (j21l) . Introducing this into the equation for 
Iq we get 



Iq = Iq{V = -SST, ST) = -nST . 



(22) 



The last equality is our definition of the thermal conduc- 
tance K. In this section we measure all energies in units 
of \U\. The dependence of k on a; = 2£d + U \s shown 
in Fig. ^ for the temperature T — 0.005 |C/| and rela- 
tively small asymmetry of the junction, parametrised by 
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FIG. 2: (Color online) Dependence of the thermal conduc- 
tance hi normalized to Hiu ~ kB\U\/h on x — {2ed + U)/\U\ 
for temperature T = 0.005![/|, ST = 0.001|[/|. Figure a) is for 
strong junction asymmetry, while b) is for very large asym- 
metry. 



the value of in the units of \U\, with F^, • F^ = 1 in 
the same units. 

Calculating the currents we have assumed symmetric 
voltage and temperature difference with Vl.r = i.V/2 
and Tl/r T ± ST/2. For positive ST the heat will nor- 
mally flow from the left to the right lead. Thermal rectifi- 
cation can be defined^Si as the dependence of the heat cur- 
rent \Iq \ or the nonlinear thermal conductance k defined 
in Eq. on the sign of temperature difference. Figure 
^ shows the effect of heat rectification in the asymmet- 
ric molecular junction with Tl/Tr = 4, the temperature 
T = 0.005|f7| and two values of ST = ±0.003|C/|. Due to 
phase space restrictions, which make the transition rates 
for two - electron processes V and ST dependent, the ex- 
cess heat current SIq = \Iq{+ST)\ — \Iq{—ST)\ depends 
on X and changes sign around |a;| w 0.022|C/| for a given 
set of parameters. 

Both the heat flux and thermal conductance are the 
symmetric functions of x for small values of asymmetry 
(see Fig. ([2^)). However, for very strongly asymmetric 
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junctions, the conductance starts to be the asymmetric 
function of the detuning parameter. This is shown in the 
Fig. (12b)- The rectification coefficient defined as 5Iq /Iq, 
changes sign as the function of x. For the parameters in 
figure (d^) it takes the maximal value of about 5%. 
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FIG. 3: (Color online) Dependence of the Wiedemann-Franz 



ratio L normalized to its Fermi liquid value Lo 



on the 



parameter x — 2ed + U in units of \U\ for T = 0.665|f/| and 
ST = O.OOllC/l. Couplings to the leads measured in units of 
\U\ are normalized by Fl • Fr — 1. 



V. WIEDEMANN-FRANZ RATIO AND 
QUANTIZATION OF ELECTRON THERMAL 
CONDUCTANCE 

It has been predicted^ that at low temperature the 
phononic thermal conductance of one dimensional dielec- 
tric wire is universally given by 



T 



(23) 



leading to the universal, material and temperature inde- 
pendent ratio kq/T = which sets a fundamental 

quantum limit on heat flow^^ in an analogy to quantized 
charge conductance. The quantized value of the heat 
conductance is expectedS^ in one-dimensional structures 
independently of the statistics of heat carriers. So this is 
valid for phonons, electrons and also particles with frac- 
tional statistics^"*. This result has been experimentally^^ 
confirmed for phonons, electrons and photons. 

It is an easy exercise to show that at T ^ the thermal 
conductance through our negative U molecule reduces to 



k{T ^ 0) 



'\u\\u\ 3 h - (1 + N) 



(24) 



Thermal conductance, which in the present context con- 
sists of only electron contribution is linear in temperature 



at low temperatures. The heat quantum, however, is not 
universal and depends on the coupling amplitudes F^/^ 
and the dimensionless distance x from the degeneracy 
point 2e(i + U = 0, measured in units of It is impor- 
tant to notice that the linear charge conductance at very 
low temperatures takes on the T independent value 

A ^2 



Go(r = o) 



Tr r, 



\U\\U\ (1 + 1x1)2 h 



(25) 



which ensures the validity of the Wiedemann-Franz (WF) 
ratio in this limit 



k{T 0) 



3e2 



— Lq . 



(26) 



G(T 0)T 

being one of the signatures of the Fermi liquid. 

For an arbitrary temperature and in the non-linear 
regime the above ratio of thermal to charge conductance 
takes on T, ST and x dependent values L{x,T) = 
The function L{x, T) in units of the Lorenz number Lq is 
shown in figure ^ iov T ^ 0.005|;7 and 6T = 0.001|C/|. 
At low X it takes values close to 1, but for larger x 
one observes strong departures from the Fermi liquid 
value Lo, like in the quantum dots with a large Coulomb 
interaction^. Its dependence on the coupling asymmetry 
is rather weak. 
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FIG. 4: (Color online) Dependence of the nonlinear dimen- 
sionless figure of merit ZT on the parameter x — {2ed + 
U)/\U\. The temperature T = 0.005|!7| and ST = -0.003\U\ 
and the couplings are normalized by Fl ■ F_r = 1 in units of 
\U\. 

We have also calculated the dimensionless figure of 
merit which indicates heat to voltage conversion effi- 
ciency of the deviceSI. It is usually defined as ZT = 
GS^T/k. For most devices thermal conductance contains 
both electron and phonon contributions. In the present 
geometry phonons do not contribute. The dependence 
of ZT on X for a few values of asymmetry is shown in 
Fig. Its value increases with the asymmetry of the 
coupling and is a non-monotonous function of x. The 
maximum increases with the anisotropy of coupling as it 
is evident in the figure. 
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VI. SUMMARY 

We have calculated the linear thermoelectric transport 
coefBcients of the single electron molecular transistor in 
the limit of large electron-phonon coupling leading to a 
negative effective charging energy. In this limit the pair 
tunneling processes contribute mainly to charge and heat 
transport. Pair dominated transport shows strong de- 
partures from the Fermi liquid characteristics. We have 
found strong deviations from the Wiedemann- Franz law 
and a non-universal value of thermal conductance quan- 
tum. The thermal quantum i.e. the ratio k/T evaluated 



at a low temperature and the Lorentz function L{x,T) 
depend on x and asymmetry of the couplings. The ther- 
moelectric figure of merit ZT characterizes efficiency of 
the device for thermoelectric applications. This parame- 
ter takes on fairly low, albeit in general x and T depen- 
dent values for the negative U quantum dot studied here. 
The maximal value of ZT is about 0.25 for the strongly 
asymmetric junction. 
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